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Abstract 

Sums of walks for charged particles {e.g. Hofstadter electrons) on a square lattice in 
the presence of a magnetic field are evaluated. Returning loops are systematically added 
to directed paths to obtain the unrestricted propagators. Expressions are obtained for 
special values of the magnetic flux-per-plaquette commensurate with the flux quantum. 
For commensurate and incommensurate values of the flux, the addition of small returning 
loops does not affect the general features found earlier for directed paths. Lattice Green's 
functions are also obtained for staggered flux conflgurations encountered in models of 
high-Tc superconductors. 
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I. Introduction 

The problem of electrons confined to a lattice in the presence of a transverse magnetic 
field is a classic problem of mathematical physics with many applications in condensed- 
matter physics. The early works of Hofstadter [1], Wannier [2], and Azbel [3] focused 
on the spectral properties of non-interacting electrons as a function of the energy and the 
parameter a = (j)/(j)o where cj) is the magnetic flux-per-plaquette and 0o = hc/e is the flux 
quantum. For commensurate a = p/q the spectrum has q subbands with special scaling 
properties as q grows larger, eventually turning into a Cantor set for irrational a. Hence, 
the behavior of the spectrum and the wave function is very sensitive to the exact value of 
a and the energy E; (the intricate structure of the bands as a function of a and E has 
been dubbed the "butterfly"). 

Lately the renewed interest in interacting electrons in two dimensions, in the context 
of the Quantum Hall Effect and various theories for high-T^ superconducting materials for 
instance, has led to more recent investigations of this general problem. Over time, these 
works have yielded an extraordinary richness using a variety of techniques including path 
integrals [4], renormalization group [5], sophisticated algebraic methods [6] closely related 
to "quantum groups," and so on. 

More recently the study of non-interacting electrons in a magnetic fleld has been 
approached from a new direction — that of "localized" wave functions. [7-10] If the energy 
of the electron does not lie within a quasiband, its eigenfunction is not a bulk (extended) 
eigenstate. However, such a wave function may be localized at inhomogeneities, such as the 
edge of the lattice (ie., surface gap states) or at isolated impurities in an otherwise ordered 
lattice. These states decay exponentially as one moves away from the inhomogeneity into 
the bulk. The effect of the magnetic field on this exponential decay (and particularly 
the associated "localization length") has been studied in the so-called "directed-paths 
approximation." [7,10] 

These calculations are based on the transfer matrix approach (reviewed below) which 
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is applicable only in this approximation. While it can be shown that this approximation 
yields a good description in the extremely localized limit in the absence of any field, the 
question is more delicate when the magnetic field is present. Indeed, to the next order, 
paths with small returning loops should be added. In the absence of a magnetic field, 
this inclusion leads to a simple renormalization of the "localization length." On the other 
hand, one finds already in the directed-paths approximation that in the presence of the 
field, the combined effects of lattice periodicity and the magnetic flux yield a complex 
interference pattern with a very sensitive dependence of the localization length upon the 
flux. [10] Therefore, the effects of including even small returning loops with their enclosed 
flux are unpredictable and potentially may be very drastic. 

Including the returning loops, as explained below, becomes even more important if 
the energy to begin with is not extremely far from the band edge and the wave functions 
are consequently only moderately localized. We thus set out in this work to go beyond 
the directed-paths approximation and to calculate systematically the full Green's functions 
which include the sum over all paths. These will be given as formal expressions which 
may be expanded systematically in various ways. Notable among these is the expansion 
in terms of the path length starting with the directed (the shortest) ones; it reveals that 
the generic features found in the directed-paths limit are not drastically changed by the 
small returning loops. In this paper we restrict ourselves to the square lattice, but using 
this approach, the sums for other lattices may be calculated as well. The expression we 
derive may also be useful to analyze properties of Hofstadter electrons within the band of 
bulk eigenstates. 

Nonuniform fliix configurations have been studied recently in the context of theories 
of strongly interacting electrons relevant to high-Tc cuprates. The most important among 
them are the staggered flux conflgurations [11]. We consequently obtain here the lattice 
Green's function for the square lattice with staggered flux as well. 

The layout of this paper is as follows. In section 2, we review the model and the 
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"directed-paths-only" transfer- matrix approach. The heart of the paper is section 3, where 
we consider the addition of "backward excursions" to the transfer-matrix method. Here, 
we find recursion relations and other expressions for the various quantities involved. In 
section 4, we examine some special cases thereof. Section 5 demonstrates how the full 
Green's function can be obtained from the quantities derived in sections 3 and 4. Section 
6 considers the staggered fiux configuration, and section 7 contains a discussion of the 
results. We also include an appendix, listing some additional results. 

II. The Model and the Approach 

One can model non-interacting electrons confined to a square lattice with a 
perpendicularly-applied magnetic field by the following tight-binding hamiltonian: 

H = e^aK + V^ala.e^^-, (2.1) 

where e is the on-site energy and V measures the overlap of wave functions centered at 
neighboring sites. The energy spectrum is the same as the one obtained using the "Peierls 

— * — * 

substitution," which replaces hk with p — eA/c in the dispersion relation obtained from 
the tight-binding approximation in the non-magnetic case. The phase = — 7ji denotes 
the Aharonov-Bohm phase an electron acquires as it moves from site i to site j; physically, 
it corresponds the line integral of the vector potential from i to j. Some freedom exists in 
the choice of gauge: To characterize a constant magnetic field, one selects the phases such 
that the directional sum of the phases 7jj around each elementary plaquette is 27r0/0o 
where (f) is the magnetic flux-per-plaquette and (f)o the flux quantum (hc/e). 

Many of the physically interesting quantities are contained in or are derivable from 
the Green's function associated with eq. (2.1). One can express the (real-space) Green's 
function f, E,(f)) between sites i and / for an electron of energy as a sum over paths 
connecting the two sites: 

paths (jk) 
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where the product is over steps comprising a given path. Each path carries a weight v^, 
where v = and L is the length of the path. Since the total number of paths on a 
square lattice grows as 4^, the sum should converge provided f < |. The model has a 
phase transition from localized to extended states occurring at v = ^ . [5] 

In the strong localization regime {v <^ the electronic wave function decays quite 
rapidly as one moves away from some initial position. In this case, Q{i, /; E, cj)) is domi- 
nated by the shortest paths — those "directly" connecting the initial and final sites. [12] 
Note that on a lattice there exists in general more than one shortest path connecting two 
sites. The presence of a magnetic field gives rise to interference among various constituent 
paths. Since it is the interference phenomena that interests us and since the difference 
in the phases accumulated between two paths is proportional to the area they enclose, we 
find it convenient to study the geometry in which directed paths can enclose the largest 
possible areas. For this reason, we begin our study with paths directed along the diagonal 
of a square lattice; we will call the directed axis t and the transverse axis x. (See Figure 
1.) Summation of paths directed along a lattice axis (rather than the diagonal) is also 
possible but slightly less convenient. 

One can employ the transfer-matrix approach to sum over directed paths. [7] If the 
vector \vt-i) represents a walk or walkers at column t — 1, then the matrix T(t, 0) operating 
on \vt-i) 

\vt) = T{t,(l)) \vt-i) (2.3) 

provides all possible extensions of these walks by one forward step and weights these 
possibilities appropriately. Consecutive operation by transfer matrices hence sums over 
all possible walks 

\vt) = T{t,4>)T{t-l,4>) ...T(l,(^) l^o). (2.4) 

Summing over paths weighted by different phase factors yields the desired interference 
effects. 
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The transfer matrix connecting column t — 1 to column t is: 



T{t, 4>) 










-it(j>-K 



-it<t>'n 











-it<t>'n 























—itrfnT 














(2.5) 



where the bonds are viewed as being directed from lower to higher t. We exploit the 
gauge degree-of-freedom to furnish a transfer matrix that depends only on t (and not on 
x). This choice is called the diagonal staggered gauge [7]. As required, summing the 
phases as one proceeds clockwise around any elementary plaquette results in 27r0; we have 
set 0o = hc/e = 1 for convenience. 

Note that in addition we have applied periodic boundary conditions along the x axis. 
(However, since we approach the problem from the localization side, we may always take 
the lateral extent Xmax of the contributing paths to be smaller than the vertical size of the 
system L_\_, and so the sum over paths will not be sensitive to the boundary conditions. 
[13]) With periodic boundary conditions, the eigenvectors of all T's are plane waves: 



i2fcj 
,i3fc J 



(2.6) 



The corresponding eigenvalues are A^^^ {t, 0) = 2v cos{k± — t(j)7r). The transverse momenta 
take on values k±_ = 27rm/Lj_, where L± is the vertical size of the lattice and m = 
0, 1, L± — 1. The eigenvalues might also be called G'j^'^{t — 1, t, (/>) for they correspond 
to the Green's function for directed walks (Fourier transformed in the transverse direction) 
joining columns t — 1 and t. 

Note that just as in the continuum model, one finds plane waves in the "transverse" 
direction, leaving an effective one-dimensional problem (an harmonic oscillator in the con- 
tinuum case). The reduction to a one-dimensional problem does not depend on the 
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restriction to directed walks since the same transfer matrix also yields a "backward" step 
from column t to t — 1. Henceforth, we need only consider eigenvalues, as all matrices of 
interest are simultaneously diagonalizable by the plane-wave eigenvectors. 

In the presence of a magnetic field, the exponential decay of the electron wave function 
away from its central position is modulated by a very intricate interference pattern. In 
the strongly localized regime, many of its properties can be uncovered by investigating the 
product obtained from the directed paths only 

To-l To To 

Gf[{0,To,cP) = J] + = n^'^-^^''^) = nh^*]' (2.7a) 

t=o t=i t=i 

where for convenience, we have introduced the notation 

Ct = cos{k± - t(t)7T); (2.76) 

suppressing the (f) and k± dependences as most of the calculations considered herein involve 
a single </> and a single k_\_. G^^(0, T^, </>) supplies the Green's function for directed walks 
that begin at the origin and end in column Tq. The properties of this product have been 
the focus of a recent investigation. [10] The transition amplitude between two sites, say 
(0,0) and {xo,To), may be obtained from them through: 



GiO,0;Xo,T,) = -3-5]Gf;(O,T„;0)e-^^--°. (2.8) 

In the absence of any field, the directed Green's function is simply a product of cosines: 

Gf[{0,To,0) = [2v cos{k^)Y\ (2.9) 

When the applied flux-per-plaquette is a rational multiple of the fiux quantum (i.e. ^ = | 
with p and q relatively prime integers), the interference-induced modulation is periodic; 
in fact, the pattern on sites of a superlattice with lattice constant q times larger than 
the original lattice mimics precisely the simpler decay pattern found in the absence of 
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an applied field. [7-10] If we restrict our attention to sites on the superlattice in the 
commensurate case, then 



One can then use the following property: 



q 



Y[2vCt 



t=i 



To/q 



(2.10) 



^ I z ' - cos(^Q«_Lj, ir q IS oaa; 
* \ 2~^+^ sin(g/cx), if g is even; 

to see the relation between the commensurate case on the superlattice and the nonmagnetic 

case. Note that in the commensurate case, v~^°G^^^o(0; ^o, |) ^ 2^°/^. 

The work of Fishman, Shapir and Wang [10] addresses the behavior of this quantity 

in the incommensurate case {i.e. when q — > oo). In this case, the structure becomes 

aperiodic. For generic irrational, it has been found that ln|v~^°G''*^| increases as [In(To)]^; 

while for an algebraic irrational \n\v~'^°G'^'^'^\ increases as In(To). 

III. Allowing Backward Excursions 

The main motivation for the present work is to consider the consequences of lifting 
the "directed paths only" restriction from the previous analysis, thereby broadening its 
scope beyond the strongly localized regime. We have opted to pursue the transfer-matrix 
method; however, this choice requires reconciling that approach, so ideally suited to di- 
rected paths, with our current concern of including returning loops. We explain first how 
the formal expressions can be obtained in an elegant and minimal way; we then proceed to 
a more detailed derivation which starts from the directed paths and adds systematically 
longer and longer returning loops. The expressions gained by this latter method will form 
the basis for our conclusions on the effects of adding small loops to fully directed paths. 
With these ends in mind, let us begin our investigation of "directed paths with backward 
excursions," which we define next. 

Let Afcj^ (t, (j)) be the eigenvalue of a transfer matrix T(t, 0) which allows any amount 
of backward excursion, that is, any number of loops of any length starting from t — 1 and 



ending at t — provided only the last step bridges columns t — 1 and t. (See, for example, 
Figure 1.) The eigenvalue \kj^ (t, (/)) might also be called G^'^ (t — 1, t, <j)) for it denotes the 
Green's function linking columns t — 1 and t that admits "backward excursions." The full, 
unrestricted Green's function can then be calculated from G|-^ (s,t, 0) as will be shown in 
section V. 

In view of the fact that all of the terms comprising A^^^ (t, (j)) contain a factor 2vCt 
corresponding to the last step, it is convenient to define {t, (j)) = Xkj^ {t, </>) /2vCt with 
this last piece stripped off. Note that the walks contributing to Xkj_ (t, (j)) begin and end in 
column t — never reaching column t. Again we could use the Green's function notation 
A/cx {i-i 4>) = ^fcx' ~ 1) ^ ~ 1) 0) which is more descriptive but a bit more cumbersome, and 
so we will put off its use until section V. 

We have found that A^j^ {t, (j)) obeys the recursion relation: 



This relation can be understood in the following way: In the expansion, the first term 
on the right-hand side (the 1) indicates the option of making no backward steps. In the 



designates any amount of backward excursion beyond t — 2 which eventually returns to 
column t — 2 (including again the option of no excursion at all) and the second 2vCt-\ 
returns the walker to column t — 1. Finally this whole process might be repeated any 
number of times as indicated by the expansion. The expansion is readily summed (at 
least formally) as a geometric series. 



-.2 ^ 



1 + 2vCt-i Afcx(t-1,0) 




(3.1a) 



second term, one of the 2vCt-i^ denotes a step back to column t — 2, then the A^x (^ — 1, 0) 





Next, consecutive A's \Xk±{t — l,^), Xk^{t — 2,0),...] can be substituted into eq. 
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(3.1b) to arrive at the following continued- fraction expression: 

1 



^^kAt,^) = 7727^2 • (3-2) 



1 



etc. 

Terminating this continued fraction at the C|_s term and re-expressing A as a ratio of 
polynomials, one finds that the numerator and denominator have essentially the same form 
except that the terms containing Ct_i are absent in the numerator. More specifically, the 
procedure leads to: 

>^kAtA) = + 0{v^^^\ (3.3a) 

Dki. [t, 0; l,s) 

where 

s s 
j=r j=r k=j+2 

- 64^' E E E CljCl.Cl, + .... (3.36) 

j=r k=j+2 l=k+2 

The problem has thus been reduced to that of understanding the properties of this 
particular function Dkj^ {t, (/>; r, s). Some of its more obvious and useful properties are: 

Dk^{t,(f);r,s) = Dk^{t + n,(f);r + n,s + n), (3.4a) 
DkAtA;r,s) = DkAt,<l>;r+l,s) - 4v^C^_,DkAtA;r + 2,s), (3.46) 
DkAt,(t>;r,s) = DkAt,(t>;r,s-l) - 4v''C^_,DkAt,(t>;r,s-2). (3.4c) 

Much is encoded in eq. (3.3). However, since our motivation is to extend the analysis 
of fully directed walks to include returning loops, we present a methodical derivation which 
does just that. Along the way useful expansions for Ajt^ (t, (p) are obtained. 

Returning then to the directed paths, a first step toward including returning loops 
would permit in addition to one forward step connecting columns t — 1 and t, one backward 
step followed by two forward steps. The modified eigenvalues (t, (p) become: 

Arjt,0) = XkAt,<f>) + [Afe,(t-1,0)]' Afejt,0). (3.5) 

10 



Let us analyze the effect of this addition on the Green's function for one plaquette on the 
superlattice in the case of commensurate Hux. As in eq. (2.7a) for directed paths, the 
Green's function is a product of these eigenvalues, and this product can be separated into 
two products as follows: 
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l[2vCt 


X 







(3.6a) 



The first product is Gf^{0, q, 2) which we considered in eqs. (2.9) and (2.10); the second 
product is also readily handled; together they yield: 

(O' 5' f ) = Gf[ (O, q, f ) X (1 - a)-"^ [l - (-l)^2a'? cos{k^q) + a^'^] , (3.66) 

where 



a = 



1 + 2^2 - VI + 4^2 
2^ 



(3.6c) 



From this expression, one can see that the main effect is a renormalization of v without 
any dependence on 0. The only structural change arises from the cos(A;_Lg) term; however, 
it is of the order times smaller than the leading behavior. Next, let us consider the 
inclusion of longer backward excursions. 

As one moves further out of the strongly localized regime, one should account for paths 
with increasingly longer backward excursions, culminating in A^^^ (t, 0) which consists of 
contributions of all lengths. A^^^ {t, cj)) can be written as: 



(3.7) 



where x\^\t,(j)) designates the contribution of length (2j + 1) which therefore carries a 



weight (2v)2-'+^. The first three are: 

Al°?(t,0) = 



{2v)Ct 

{2vfCtCU 



;(2), 



(3.8a) 
(3.86) 
(3.8c) 
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We provide expressions for A^^^ {t, </>) for j = 0, 1, . . . , 5 in the appendix. They become 

cumbersome rather quickly as j increases. For j > 1, \^^^{t, (f)) can be shown to obey the 
foUowing recursion relation: 



E 



A£^(t-1,0). 



(3.9) 



Again we find it convenient to strip off the last step of the walk and study Xkj_{t, (f)). 
After generating A^^^(t, 0)'s from the recursion relation, we have observed that they could 
be expressed in the following way: 



E 



n 

i=l 



Hi + rii+i - 1 
rii+i 



rii 



(3.10a) 



where Ui = 0,1, j for i = 1,2, . . . ,j and where 



Hi + Ui+i - 1 

rii+i 



1, 
0, 

{rii + Ui+i - 1)! 



if Hi = and Ui^i = 0; 
if = and nj_|_i > 0; 

otherwise. 



(3.106) 



(n,+i)!(n,-l)!' 

The restriction on the n^'s in eq. (3.10a) ensures that the contributing walks have the 
appropriate length. Note first of all that each cosine term appears an even number of 
times, making the round trip possible. Moreover, the terms are "connected;" for instance, 
if a term contains a factor Cfs, it must also contain factors Ct-2 and Ct_i. Finally, 
the numerical component (the degeneracy) factorizes into terms which depend only on the 
number of steps in two neighboring columns {rii and rii+i). 

Expresson (3.10) may have certain practical advantages over eq. (3.3) in that it 
selects out walks of a particular length and that it is not written in the form of a quotient. 
Nevertheless, we proceed now to use it in a second derivation of eq. (3.3). Performing the 
sum of A^-'^ {t, (j)) over j 

oo 

^ (3.11a) 
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actually simplifies matters as it eliminates the restriction on n^'s, leading to 

(3.116) 



{ni,n2,.--} 



Next carry out the sums over consecutively. Using repeatedly the following sum 

A; - 1\ ^ X 



+ = l + 7T^. (3.12) 



i=o 
leads to: 



+ 



(1 - 4^2C2_^) (1 - Av^CU) (1 - ^v^CU - ^v^Cl,) 
+ 0(^8). (3.13) 



Finally finding common denominators for the first two terms, then the first three 
terms, etc., yields expressions like: 

and so on. Hence, one arrives at: 

hAtA) - '^'t'l^'l + 0{v^^^% (3.15) 

Dk^ {t,(p; l,s) 

which is identical to eq. (3.3). 
IV. Special Cases 

In this section, we consider the evaluation of Xkj^ {t, 0) in three special cases: 1) when 
V is small (which provides a straightforward extension of the directed-paths-only analysis 
in the strongly localized regime); 2) when the magnetic flux-per-plaquette is a rational 
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multiple of the flux quantum; and 3) when is small (weak magnetic field, for which one 
should be able to make some correspondence with the continuum model). 

iv.a. When v is small 

The first special case we shall treat is when v is small. Under such conditions, it is 
reasonable to view the various quantities as power series in v. We have already seen in 
the previous section that 



Al°?(t,<^) 



= 1 



(2vfcU 

(2v)^ Cf_i + 2C^_iC^_2 + C^_iCf_2 + Cf_iCf_2Cf_^ 



and so on. From these, we see that Xkj^ {t, (j)) can be re-expressed as 



(4.1a) 
(4.16) 
(4.1c) 
(4.1d) 



(4.2a) 



(4.26) 



(4.2c) 



Expressions (4.2a), (4.2b), and (4.2c) account accurately for backward excursions of length 
2, 4, and 6, respectively. 

Note that the expansion of the ln[Xkj^{t, (p)] (the argument of exp{} above) looks 
rather like the expansion of Xkj_{t,4>) itself (eq. (3.11b)) except that coefficients have an 
additional 1/ni factor 



\ii[XkAt, </>)] 



{ni,n2,...} 



(4.3) 
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and, of course, the ni = term is absent. We can verify this expression with repeated 
use of the relation 



i=i 



1 fj + k-1 



k 



X-' = 



-ln(l - a;). 



if /c = 0; 



which leads to 



(4.4) 



\n{XkAt,<l>)) = ln[i^fc^(t,(/>;2,s)] -ln[iPfc^(t,0;l,s)] + Oiv^'+^), (4.5) 

which is simply the logarithm of eq. (3.3a). 

The form given above is particularly convenient for calculating G^-^ (O,To,0), the 
Green's function for long walks that allow backward excursions, because G|-^-(0, T^, 0) 
factors into G^^(0, To, 4>) and a product of A's 

To To 

Gt"(O,To,0) = n >^kAt,<t>) = Gf;(O,T„0)xn Afejt,0), (4.6a) 
t=i t=i 

and the exponential form converts that product into a sum: 



G^-^^-(O,To,0) =G£^(O,To,0) 

X expjx; [(2^)'c,^-i + (2^)' + CUCU] + 



t=l 



(4.66) 



Let us return to the example considered in the previous section (eqs. (3.5) and (3.6)) 
in which we have examined the effect of short backward excursions (of length 2 only) 
on the interference pattern generated on the superlattice when the flux-per-plaquette is 
commensurate {(f) — |). We are now in a position to extend that analysis. Substituting 
To = q and = | into eq. (4.6b), we find 



xexp (2^)2| + (2^)4[l| + |cos(^) 



+ 



(4.7) 
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where we have used the following relations for = | 

E^'-i = 1 for ?>2, 
t=i 

Y^CU = y for g>3, 



for g > 3. 



(4.8a) 
(4.86) 
(4.8c) 



The exponential part represents the correction to the directed-paths-only approach. We 
see the same pattern here as before — that on the superlattice (Tq = q) the transverse 
momentum does not enter the corrections {i.e. there are no structural changes) at 
order t;^ or order provided g > 3, suggesting that k± enters only at order v^'^. Some of 
the sums needed to carry out this analysis to higher order are: 

q ( 2ra' 



EC?" 



^3 - 



22m+2n-l 



E 



22r 

2m 

171 + I 



m 



2n 
n-£ 



if q > m, 



cos 



1 /2m 



m 



2n 
n 



(4.9a) 



(4.96) 



ii q > m + n. This last calculation is really a mixture of two special cases — small v and 
commensurate flux {(f) — |). We will now turn our attention to the commensurate case. 

iv.b. Commensurate flux 

In the case in which the flux-per-plaquette is a rational multiple ( g) of the flux quan- 
tum, the terms in the continued fraction expression (eq. (3.2)) begin to repeat themselves 
after q terms. For instance, for = |, one has: 

1 



'^fc.L {t, I) 



(4.10a) 



1 - 



1 - 4v^Cl:,XkAt,k) 



Simplifying the fraction yields: 
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(4.106) 



Performing the same procedure for any = | , we find: 



^ Dk^ (t, g; 2, g - 1) - Av'Cl^Dk^ it, g, 2, g - 2) A^, (t, g) 

^ ^ Dk^ (t, g; 1, g - 1) - 4^2C,2_,Dfe, (t, g, 1, g - 2) A,, (t, g) ' 



(4.11) 



From here we conclude that Afc^^ {t, g) is the solution of the following quadratic equa- 



tion: 



AJ^,g) 



A..(t,g) 



+ i3,,(t,g) A,Jt,g) + C,,(t,g) = 0, 



(4.12a) 



where 



A,(t,g) 



= 4^2C2_ n,^(t ,f;l,g-2), 



i3fc,(t,g) = -i?fejt,g;l,g-l) - 4^^Ct,i?fc,(t,g;2,g-2) 
CkAt,^) =i^fc,(t,f;2,g-l). 



(4.126) 
(4.12c) 
(4.12d) 



The solution of the quadratic equation is: 



Bfe,(t,g) - ^B^,(t,g) - 4 AJt, g) Cfe Jt, g) 
2 AJt,f) ' 



(4.13) 



We have used this last expression to calculate Xkj_{t, g) for relatively small g. Some 
algebra and trigonometry reveal that: 



Afe^(t,l) = 



1 - 4^2 + 8v^Cf_2 - ^J{l-4v^)^ -16v'^sm^{2k±) 



1 - + 8v^C^_3 - 7(1 - 6^2)2 - 16v^cos'^{3ki_) 



(4.14a) 

(4.146) 
(4.14c) 



One can see a pattern arising from these expressions. For instance, we can see here and 
show in general that the discriminant (B^ — AAC) is independent of t. 
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Note that q = 1 corresponds to one flux quantum per plaquette or equivalently no 
field at aU. The smaU field limit, on the other hand, would require taking q to infinity. 
Fortunately, this limit is accessible by other means. 

iv.c. The small field limit 

The third special case we consider here is when the magnetic field is small. Toward 
this end, let us perform a Taylor expansion of Xk±{t, (f)) around {k_i_ — t(f)7r): 



(4.15) 



i=0 



It is our goal in this section to calculate the z = 0, 1 and 2 terms A[.*^ (t, (/>) above 



Instead of expanding A^j^ {t, (f)) directly, we will begin by expanding A^^^ {t, 0). Recall 



ing eq. (3.10a) for A^^^(t, 0), we clearly need an expansion of C|_„- around {k± — tcpir): 



= CI - Mi^ct's, + 



i)ct'sl 



ecf] + 



(4.16) 



where St — sm{k± — t(f)iT). Below we present the expansions of A^"^^(t, 0) to order 0^ for 
J = 0,l,...,6: 



2 



={2vy Ct - 27T<l)CtSt + -i-{2St - 2Ct) 



^{2vf 2Cf - 107T(j)CfSt + 
=i2vf hcf - U7T(l)C^St + 



n 



2m2 



2 



--{2vf 14Cf - 1867r0Cj5t + ^^-^{22AQC^Sl - 374Cf ) - 



2^.2 



= i2vr 



2^2 



42C/'J - 7727T(l)C^St + -^{13348C°Si - 1748C/") - 
132Cf - 3172TT(t)C}^St 



+ ^^(73340Cf <S| - 7916Cf ) - ... 



2 J,2 



(4.17a) 
(4.176) 

(4.17c) 

(4.17d) 

(4.17e) 

(4.17/) 

(4.17^) 
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Next, we collect terms of the same order in and resum those series. Toward that 
end, we define: 

oo 



Ct 



.(2) _ ^ 



From eqs. (4.17a-g), we can conclude, for instance, that 

[Pj] = { 0, 1, 5, 22, 93, 386, 1586, ... }. 

These four series are summed as follows: 

1 



2x 
1 



j=0 



3-^ - 2 



l-(l-4a;)V2], 
(l-4a;)-i -(l-4a;)-V2j^ 



^g^. = -(l-4x)-5/2-2(l-4x)-2-2(l-4x)- 
+ (l-4x)-i + J(l-4a;)-i/2, 



-3/2 



j=0 



■3/2 



Using these we arrive at: 



1 - yi - iQv^c^ 

8^;2^2 : 



and 



(4.18a) 
(4.186) 
(4.18c) 



(4.19) 



(4.20a) 
(4.206) 



(4.20c) 
(4.20(i) 



(4.21a) 
(4.216) 



2„.2r7-— 3 J 0/02 /^2^ 



5 + C/ 



1-C/^ 



(4.21c) 
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where U = {1 - IGv'^Cf)-^/^. 

Note that the hnear term breaks the (p ~0 symmetry. Two comments should be 
made: (i) (f) —(f) with x — > —x (or equivalently k± —k±) is stiU a symmetry and (ii) 
most importantly all physical properties which depend on the absolute value \G\ will be 
invariant under — > —0 or a; — > —x separately. This is consistent with the expectation 
that no Lorentz force (and no parity symmetry breaking) should occur for states which 
do not carry current density. Such effects will be seen only if the dependence of G in the 
t direction goes like e*'^!!* so that the states have a nonvanishing current density in the t 
direction. [14] 

V. The Full Green's Function 

In this section, we show how to obtain the full Green's function Qk^ (0, To, (j)), which is 
the appropriately weighted sum of all walks connecting columns and Tg. (Recall that the 
2D problem was reduced to ID.) Thus far, we have calculated X^^ {t, 0) = G'j^^- (t — 1, t, 0), 
which includes walks which advance from column t — 1 to column t with any amount of 
backward excursion prior to the forward step. In this section we will find it more convenient 
to use the Gfe^ (t — 1, 0) notation. We have also already seen the product: 

To — I 

Gt^-(0,To,</>) = n Gt^-(t,t + l;</>); (5.1) 

which supplies the Green's function for the restricted walks that begin at the origin and 
end in the column, again with the last step being the only one connecting the To — 1 
column to the column. This is the Green's function for the first time the walker lands 
in column Tq. 

Now we must allow the walker to go beyond Tg or to turn back and return later to T^. 
Let us denote this sum by Q^,'^^'' {To,To; </>); it is the Green's function for all walks starting 
in column Tg and returning to that column. Note that ^^"'^(Tq, Tq; (j)) obeys the following 
recursion relation: 
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glf{To,T,;cl>) = 1 + 4v^Cl Gl'^iTo - 1,T, - 1,0) gif{To,ToA) 

+ 4^'C|^+i G{;^"(T, + 1,T, + 10) gif(To,To;<t>), (5.2a) 

where G-^'^' is the same Green's function as we have aheady found except that the initial 
position is to the right of the endpoint. Therefore, the formal solution is: 



6?,^^"(T„ T,; 0) = yi-Av'C-'^G'^^^iTo - 1, T, - 1, 0) -4^^C^^+iG^f (T, + 1, + 1; 0) 

(5.26) 

Finally, the full Green's function ^^^'^(0, Tq; 0) is simply the product: 

gif{Q,T,-A) = Gl^^{Q,To-A)Qlf{To,To-A), (5.3a) 

To — 1 

n Gt^-(t,t + l;0) 

^{^"(0, T„ 0) ^° 



1 - Av^ClG'.^iT, - 1, T„ - 1; 0) - Av^Cl^^Gi^- {T^ + 1, T, + 1; 0) 

(5.36) 

This last equation is the central equation of this paper: It gives the formal expression 
for the Green's function of an electron on a lattice in a magnetic field. The dependence on 
X can be obtained by an inverse Fourier transform with respect to exp{ik±x). In principle, 
all physical properties may be obtained from this expression. This is certainly true for 
E outside the spectrum. For the energies within the "butterfly," one may still use this 
formula for walks which do not reach the boundaries of the system (or wrap around the 
torus for if periodic boundary conditions are assumed, other delicate questions then arise 
[13] if the two sizes L± and are not integer multiples oi q.) 

In the next section, we look at the Green's function for charged particles propagating 
within the staggered-fliix configurations that have been studied in the context of high-T^ 
superconductivity. 

VI. Staggered Flux Configurations 

21 



Beginning with the work of Affleck and Marston [11], the notion of staggered flux 
configurations was introduced and studied in the context of theories of high-Tg supercon- 
ductivity. It may, therefore, be useful to have expressions for electron propagation in such 
flux configurations. On a square lattice, such a configuration will be given by a chessboard 
arrangement of interpenetrating sets of plaquettes with (say) +(f) on the white squares and 
—(j) on the black ones. (See Figure 2.) The most studied case is = tt, for which the 
phase factors acquired around the plaquettes are ±1, so that time reversal symmetry is 
unbroken. The sum of walks in this particular case was obtained by Khveshchenko, Kogan 
and Nechaev [15]. In the present section, we obtain these sums for arbitrary staggered 
flux ±0. 

First, we must choose a gauge: this time we assign the phases only to the vertical 
bonds of a square lattice. Considering them as directed upwards, we assign a factor 
7 = e*"^/^ for bonds pointing from sublattice A to sublattice B (shown as dashed in Fig. 
2) and 7* = to the bonds pointing from sublattice B to A (solid in Fig. 2). So a 

staggered arrangement of phases ±0/2 is assigned to the (upward-pointing) vertical bonds. 

Next, we assign a factor a {ol~^) to a step in the positive (negative) x direction 
and similarly (3 (/3~^) for the y direction. A natural choice is a = e*'^^ and (3 = e^^y 
which yields the Fourier transform of the generating functions. We define the following 
generating functions: first, the local generating functions for a step 

A^B: QA ^ a + a"^ +7(/3 + /3-^) 

B^A: QB = a + a-^ + -f-\p + p-^), (6.1) 

one for each sublattice; then a global generating function for all paths of length 2L: 

G^^^^{(f)) = g\gB = 4:{cos'^ ka; + cos'^ ky + 2 cos cos ka; cos ky) . (6.2) 

The probability amplitude to reach a given site (x, y) is the coefficient of o;^/?^ (up to a 
trivial normalization) or the inverse Fourier transform: 

Vt\^) « \[ \[ dk. dky Gl'l\ct^) e-'^^-^^yy, (6.3) 
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where the system size is assumed to be infinite and we have thus gone over to integrals 
over the momenta. Note that in order to obtain the normahzed Green's function, qa and 
gs must each be multiphed by | and G2L by 2~^^. 

Of particular interest are the expressions for walks which return to their initial point. 
The weighted closed magnetic walks are given by the coefficient of q;°/9° = 1 in the product: 



(6.4) 



where 7 = e*"^/^. Toward this end, we expand the above products as follows: 



(A) 



(6.5a) 



(a + a-i) + 7-^(/3 + /3-^) " = J] ( f^^ V« + c"^)^^"^ [7-^(/3 + Z?"^)] ''^ (6.56) 



-(B) 



where l'^^^ is the number of horizontal (vertical) steps initiated from the A sublattice 

and likewise for the B sublattice. Consequently, L^f^ = L^x^^ — L. Next, we 

expand the remaining a products, introducing the variables and which represent 
the number of steps in the positive x direction initiated on the A and B sublattices, 
respectively. The coefficients of a will thus be: 



E ( ,(_Aj +^ ^ 



t{B) , 

JB) l« - 



(6.6) 



To obtain the q;° term, the following relation must hold: 



(6.7) 



where is, of course, the total number of horizontal steps. Note that must be even. 



Replacing by ^ — , we find that the coefficient of is 







4"' ,^(.4) 



E 



X 

{A) 

(A) ^ '+X 
+x 



J^x 



(A) 
+x 



(6.8) 
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Repeating the same procedure for the y direction and combining the results yields: 



where we have used L^"^^ = L — L^^-* and Ly = 2L — L^- 

Using eq. (6.9) we have found G^*^^ for several L's. They are: 

G^*"^ = 1 (6.10a) 

G^*"^ =4 (6.106) 

^s^tag = 28 + 8 cos (6.10c) 

G'^*"^ = 256 + 144cos(/) (G.lOd) 

Qstag ^ 2716 + 2112 COS + 72 cos 20 (6.10e) 

Gfo"*^ = 31504 + 29600 cos + 2400 cos 20 (6.10/) 

Gll""^ = 387136 + 411840 cos + 54000 cos 20 + 800 cos 30 (6.10^) 

Cf/s = 4951552 + 5752992 cos + 1034880 cos 20 + 39200 cos 30 (6.10/i) 
Glf^ = 65218204 + 80950016 cos + 18267200 cos 20 + 1191680 cos 30 

+ 9800 cos 40 (6.10i) 
Glf^ = 878536624 + 1148084928 cos + 307577088 cos 20 + 29070720 cos 30 

+ 635040 cos 40 (6.10j) 
Gfo"*^ = 12046924528 + 16407496800 cos + 5030575200 cos 20 + 625312800 cos 30 

+ 24343200 cos 40+ 127008 cos 50 (6.10A;) 
GgL*^ ^^^^ ^® obtained as the coefficient of v'^^ in the following expression: 

— J dk±^l-8v^ - 8v^ cos(0/2) cos(2A;i) + 16w^ sin2(0/2) sm'^{2k±_) | ' (6.10£) 
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As we work toward procuring a more general expression, let us concentrate on the 
following sum: 



(6.11a) 



which is the latter portion of eq. (6.9c). Introducing L^^^ = Lr^ — L^^^ this may be written 



as: 



L L 



t{A) j{B) 



X X 



1 
1 

2^ 



7 



t(A) ,(S) 
^x 



27r 



de 



1 + 76 



ie 



1 + 7-ie'^ 



— / dOe'^^^-^-^ 2 cos ^+ (7 + 7-^) 
2tv ,/o L -I 



27r 



COS ^ + cos 



2 J 



(6.116) 



We may now use the identity: 



1 

2^ 



27r 



de 



11 + — 1 cos 9 



" / cos m0 



(6.12) 



^sin m9 J (n + m)!\ 
We identify cos ^ = /// -yZ/x^ — 1 which implies fi = ±i cot ^ such that the integral becomes: 

2^ 



27r (^2 _ 1) 
Therefore for the sum 

Lx 



1 /"^"^ r 1 

y e*^^-^-^-) [yu + Vl^^ - 1 cos 



(6.13) 



E 



(A) 



{A) 7 



2' (2L-Lx)\ ^ 



The final result is: 



Lx 
2 



Pi^"^^^(icot^) 



(2L-L,)! ^ 



(6.14) 



(6.15a) 
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with 



C{L) = (-i)*(sin^) ^2^L! (6.156) 



We find agreement with the known cases [15] of (/> = for which eq. (6.14) yields (^^) 
and = TT for which it gives ( ) . 

One can generahze the expressions to walks that do not close. To reach a given point 
{xo.Uo) the only diff'erence will be in the summation over l^J since what we need now is 
the coefficient of or 

201^ + 4?) -(4^^ + 4^^) = -o. (6.16) 
Therefore l[^J = ^ + ^ - l^_^J and the sum is: 

/4^N/ 4"^^ \ f \ 

,(A) V'+a; / \ 2 2 / V 2 ^ 2 / 

Likewise (to get the coefficient of /?^°) the summation in the y direction yields: 



2 2 



and the final expression changes to: 



2L ^ 



GT = C(L) E (./,%. K . ) T^TF^^-f -'-'(icot --). (6.19) 



L^=0 V 2 ' 2 / \^ 2 ' 2 

Alternatively, one could choose a = e*'^" and (3 = e^^y and study G2L(a, /3) which can then 
be Fourier transformed back into real space. 

VII. Discussion 

In this paper we have derived expressions for the square-lattice Green's function of a 
charged particle in a magnetic field. We have also obtained similar expressions for the 
so-called staggered flux configurations with arbitrary flux ±0. The results were obtained 
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within the Peierls ansatz in which the magnetic field is represented by an extra phase 
acquired by the nearest-neighbor hopping term. 

Our starting point was the directed-paths-only approach initiated in the study of 
strongly localized wave functions. From there, we have shown how to include system- 
atically the returning loops and obtain the Green's function for walks with backward 
excursions, and then we have demonstrated how the full Green's functions were realized 
in terms of the former. 

Our main goal was to see if the inclusion of small returning loops would modify in 
some essential way the behavior found in the directed-paths approach. The analysis of 
eqs. (3.5), (3.6) and (4.1) through (4.8) shows that the answer is negative in the sense 
that the product over cosines with changing phases obtained for directed paths is modified 
into a product of more complicated trigonometric expressions containing the same cosines. 
Of course, the behavior will change but the essential features: strong sensitivity to the 
value of the field, simple scaling in q for the commensurate case 4> — p/q, and an aperiodic 
but deterministic behavior for irrational (f), all are going to survive the addition of the 
small returning loops. Repeating the analysis done in ref. 10. for the product of the new 
trigonometric expressions would be a formidable challenge, since even for the former case, 
the analysis was a substantial and sophisticated task. 

Our expressions might aid in the study of other interesting issues as well: For example, 
the transmission through a slab of finite width in which the wave function does not decay 
all the way to zero, might also be studied. An especially interesting question concerns the 
emergence of a Lorentz force (and the associated parity symmetry breaking) in the case 
which deals with such a current carrying state. 

Another interesting issue is how the continuum limit may be approached: Naively, 
one would expect that if the magnetic flux is small such that the magnetic length is much 
longer than the lattice spacing (but smaller than the system size) the lattice results will 
coincide with those obtained in the continuum. In particular, the exponential decay should 
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change into a Gaussian ~ exp{— c5r^} behavior. [4] 

Finally, these Green's functions should be correct (analytically) for energies within 
the band. One should, of course, add a small imaginary part to the energy (and de- 
fine advanced and retarded Green's functions) to deal with the singular behavior at the 
eigenenergies. In addition, the question of boundary conditions for walks hitting the 
boundaries should be handled with special care. In principle at least, all physical infor- 
mation could be extracted from these Green's functions. It will certainly not be straight- 
forward to extract physically relevant information from the formal expressions, but in view 
of the importance of such information for many physical problems of interest (especially 
the many-body extensions studied in the context of the Fractional Quantum Hall Effect, 
high-Tc superconductivity, anyonic physics, etc.), it may be worthwhile to pursue those 
directions. 
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Appendix 

In section 3 we have introduced the function Xkj^{t, 4>), the eigenvalue of the transfer 
matrix T{t,4>) which aUows any amount of "backward excursion" foUowed by a single 
forward step. We called ^^^{t, 0) the contribution with path length {2j + 1). Using the 



recursion relation (3.9), we have generated A^''^ {t, (f)) for j = 0, 1, . . . , 5. They are: 



{2v)Ct 
{2vfCtCU 

/o6 , o/^4 /o2 , /o4 , /o2 /o2 /^2 

C8 1 /^2 1 1 o/^4 /^2 /^2 

, /o2 /o6 , o/^2 /o4 /o2 , /o2 /o2 /o4 



= {2vyct 



{Al.la) 
{Al.lh) 
{AI.lc) 
{Al.ld) 



't-3 
3 



, /^2 y^2 /^2 /^2 
-I- *^t-l<-^t-2'-^t-3'-^t-4 J 

+ ^C^_iC^_2 + QC^_iC^_2C^_^ + 2C^_iCf_2Cf_ 

I Q^4 ^2 ^2 ^2 , ^2 ^8 , Qr'S ^2 
+ ^'-^t-l"-^t-2"-^t-3"-^t-4 + *-^t-l*-^t-2 + ^*-^t-l*-^t-2*-^t- 

I Q/^2 /^4 /^4 1 0/^2 /^4 /^2 /^2 1 /^2 /^2 
+ o<-^t_i<-^t_2<-^t-3 + ^<-^t-l<-^t-2<-'t-3'-'t-4 + '-^t-l'^t- 

I o/^2 /^2 , ^2 ^2 ^2 ^4 

+ ^«-^t_l«-^t_2'-^t-3'-^t-4 + ^t-V-^t-^^t-Z^t-l 

I /^2 /^2 /^2 /^2 /^2 

where Ci = cos^A;_l — tcfrn^. These expressions can be seen to obey eq. (3.10a). 



(Al.le) 



y^6 



(Al.lf) 
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Figure Captions 



Figure 1. The diagonal lattice with t and x axes labeled and an example of a directed 
path with backward excursion. 

Figure 2. The Staggered Flux Configuration. 
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